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In this paper we consider dark matter particle annihilation in the gravitational field of 
black holes. We obtain exact distribution function of the infalling dark matter particles, 
and compute the resulting flux and spectra of gamma rays coming from the objects. It 
is shown that the dark matter density significantly increases near a black hole. Particle 
collision energy becomes very high affecting relative cross-sections of various annihilation 
channels. We also discuss possible experimental consequences of these effects. 

Keywords: Dark Matter; Black Hole; Gamma-Ray Astronomy. 



1. Introduction. 

The Dark Matter abundance of the Universe is now estimated as ~ 23% of the total 
energy contenfM The most popular hypothesis about its nature is that it consists 
of heavy stable elementary particles surviving from the early Universe. There are 
strong cosmological reasons to believe that the particles are very cold now (« < c). 

In spite of no lack of theoretical candidates for the role of dark matter particle 
(hereafter DMP), its nature is still hazy. A very promising way to make the situation 
clearer is to observe astrophysical annihilation of the DMPs. In this article we 
consider dark matter annihilation in the gravitational field of a black hole. There 
are several reasons why this task is important. First of all, the greatest dark matter 
density in our Galaxy should occur in the center, and then the annihilation rate there 
must be the highest. A possible signal from the central source is a subject of wide 
speculation^. Experimental data indicate a supermassive black hole in the centre of 
the GalaxjP. Secondly, the dark matter density should significantly increase near 
a compact massive object (see below). Thirdly, as it will be shown in this article, 
some very specific physical effects can appear nearby a black hole. 

The influence of the black hole gravitational field on the dark matter annihila- 
tion has already been considered-^, but the physical models were rather rough. In 
particular, they completely disregarded a significant increase of the DMP velocities 
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and strong anisotropy of their phase space distribution. Meanwhile, as it will be 
shown below, the accreting DMP velocities are comparable with the speed of light. 
If we surmise the DMP mass to be equal 40 GeV, the particle collision energy can 
reach hundreds of GeVs. At such high energies the cross-sections of various annihi- 
lation channels vary significantly. This can strongly affect the annihilation rate and 
annihilation products. 

The plan of our work is the following: first of all we obtain the exact phase dis- 
tribution function for noninteracting particles accreting on a Schwarzschild (nonro- 
tating) black hole in the spherically symmetric case. We assumed for simplicity that 
the particles far from the black hole are nonrelativistic (v -C c). This assumption 
is well motivated in physically important cases, as it will be shown below. Then we 
calculate the number of collisions and the product distribution, and finally compute 
the product outlet to the distant observer. In conclusion, we discuss the obtained 
results. 



2. Dark Matter Particle Phase Distribution in the Vicinity of a 
Black Hole 

The gravitational field of a non-rotating black hole of mass Mbh in vacuum can be 
described by the Schwarzschild metric: 



dr 2 
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where r„ = - and c is the speed of light. The trajectory of a particle falling 

c l 

onto a Schwarzschild black hole can be written (Ref. [51 equations 101.4 and 101.5) 
as: 
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where Eq, m and M are the total energy, mass and angular momentum of the par- 
ticle respectively; E and M are, of course, integrals of motion. As we have already 
mentioned, there are strong reasons to believe that the dark matter particles are 
cold. Then their speed is mainly determined in the Galaxy by its gravitational po- 
tential and can be estimated as a few hundred kilometers per second. Consequently, 
the considered particles are nonrelativistic far from the black hole which allows us 
to simplify the equations. Besides, we introduce here the system of units in which 
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the speed of light and the gravitational radius of the black hole are equal to the 
dimensionless unity: c = 1, r g — 1. Then equations (12IM[) can be rewritten as: 



f = / ^^=^= ( 4 ) 



° <lr (5) 



Ff v/ 'r 2 — a 2 (r — 1) 

where a = — . For the radial component v r , tangential component v tan and the 
module v of the particle velocity @ we obtain: 



\Jr 2 - a 2 (r - 1) a . . ... 

Vr = 7= ; Vtan = ~\ (6) 

ryr r V r 

Hence, the angle 6* between the direction to the center of the black hole and the 
particle trajectory is: 



r r 
From ([7]) one can see that all the particles reaching a radius r have the same velocity 
module. It is convenient to introduce the quantity N of particles in a unit volume 
moving at the angle 9 to the direction to the center of the black hole in a unit solid 
angle. Then the quantity of the particles crossing the sphere of the radius r in a 
unit time in a unit solid angle in the phase space i:Q 

4nr 2 H v cos 9 drdVl (9) 

If r is fixed, we obtain from ([8]): 

cos6>dft = cos(9d(27r(l-cos60) =7rd(sin 2 (9) = Trrf " ^ = i^—^- da 2 (10) 
Besides, we substitute the interval of physical time dr for the interval of the universal 



Schwarzschild time coordinate dt — drJ 1 — £ . Then we can rewrite (j9|) as: 



47r 2 N (viZlLdaadt (ii) 
r 

Let us consider the particles from the interval [a; a + da\; a is an integral of motion, 
so a and da remain constant along the bundle of trajectories. Since the system is 
stationary, the quantity of the particles considered crossing any sphere of the radius 



a By v r , vtan, and V we imply the true, physical components of the velocity, in contrast to the 

y/—g r r dr dr 

coordinate ones. For example, v r is equal to instead of — 

xFglt dt dt 

b We use that in the Schwarzschild metrics the surface of of a sphere of radius r is equal to 47rr 2 . 
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r (if they cross it at all) in a unit of the Schwarzschild time t must be equajfj. Fix 
(fTTjl we obtain: 



4tt 2 H 1)3 = const (12) 

r 

In order to find N, we should impose a boundary condition on the infinity. From 
(fT2| it follows that H — > if r — > oo, while it would be reasonable to expect that 
N — > rioo/Air, where rioo is the dark matter particle concentration at infinity. The 
contradiction is a consequence of our assumption that the velocity of the particles 
at infinity Voc is zero. According to Ref. [5l the gravitational capture cross-section 
of nonrelativistic particles is: 

2 



4-rrr 2 . 



c 



The cross-section becomes infinite if — > 0. In order to avoid this difficulty, 
we should consider that in practical situations is never equal to exactly: for 
example, near the Earth the dark matter particle speed is Woo ~ 300 km/s. We cut 
the solution (fi"2|) on the radius r m , where the velocity given by (|7|) becomes equal to 

Woo. According to l|7|). r w = —5-. The dark matter distribution at larger distances 

will be considered undisturbed by the gravitational field of the black holcE). Then 
the quantity of the particles crossing the sphere of the radius in a unit time in 
a unit solid angle in the phase space is: 

47rr^ c iiooVoo cos - — dr = n^Voo cos d£l dr (13) 

47T 

t — 1 ndcP" 
According to (fTO]) cos 0dfl = tt — ^— da 2 . As 7*00 ^ 1, we obtain cos 6 dfl = ; 



dr = dt. Then we can rewrite (|13)) as: 

^ala 2 dt (14) 

Woo 

Comparing equations (fl"2"|) , (113(1 and (|14p . we derive that the constant in (fT"2"l) is 
equal to: 



COnst = 7T- 



Substituting it into (I12p . we obtain 



(v/-i) 3 (15) 



c We use also that the Schwarzschild metrics is static. 

d This supposition is true if Uoo c. In fact, if Voo = 300 km/s, then roo = 10 6 r 9 i.e. gravitational 
field on the radius roo is weak 
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Deducing the initial formula (|12|) and considering the particles crossing the sphere 
r at the angle 8 to the direction to the center, we imply that their world lines begin 
at infinity. It is true for the angles 



cos > 4 ( r R S ( 2 - r) (16) 



where 



x > 

E(x) = <( 0, x = 
-1, x<0 

(see the details about particle trajectories in the Schwarzschild field in Ref. [5]). The 
world lines, for which the condition (fl7j|) fails, should begin on the event horizon. 
Of course, there are no particles in these angles. Finally, we obtain for H: 



Jr 2 -Mr- 1) 
cos6» > i 5(2 - r) 



47TUoo (y/r - l) 3 ' r ' ny\ 

0, cos0<^- -i- =iS(2-r) 

r 

This is the formula which describes the collapsing particle distribution around the 
black hole. 



3. Calculations and Results 

Using the distribution function (TTT)) of DMPs, one can calculate the annihilation 
signal (measuring by a distant observer) for various annihilation channels and var- 
ious dependencies <j(E). As an illustration, we consider the following important 
model process: a particle and an antiparticle of the dark matter annihilate into 
two photons(xx ~~ * ^l)- The DMP mass is assumed to be todm = 40 GeV. The 
cross-section is described by the Breit-Wigner formula: 

w^www (18) 

Here E is energy in the cm. system, E r is the resonance energy. We take E r = 
91 GeV, r = 2.5 GeV which approximately corresponds to the values of Zq boson, 
though we do not propose any concrete annihilation process, considering just a 
model problem. We suppose uo = 10~ 4 GcV~ 2 in order to provide a characteristic 
weak cross-section out of the resonance er(50GeV) ~ 10 -8 GeV -2 . The photons 
were assumed to outcome spherically symmetrically in the cm. system. The photon 
energy, of course, increases with the collision energy. Besides, the centre of mass 
of the colliding particles is generally speaking not at rest. So in the laboratory 
reference system the photons are created asymmetrically, and their energy depends 
on the angle of producing. 
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We use stationarity of the problem to calculate the photon outgoing to infinity. 
The general theory of relativity effects are taken into proper account. First, the 
photons undergo redshift: a photon produced on the radius r with the energy Eq 
outcomes to a distant observer with energy 



are captured by the black hole i 

We assume the black hole mass to be equal Mbh = 3T0 6 Af© (which corresponds 
to r g = 10 12 cm), v x = 300 km/s, and = 2 • 10~ 2 cur 3 (~ 0.8 GeV/cm 3 ). 
The last value is characteristic at the Sun system vicinity, for the Galaxy centre 
it must be underestimated. Unfortunately, there is no reliable estimation for the 
dark matter density in the centre of the Galaxy. The results of the calculations are 
represented in Fig. Q] 




(19) 



Second, the photons produced at the angles 
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Fig. 1. The spectrum of gamma-ray emission produced by the dark matter annihilation into two 
photons. The cross-section is discribed by | |18| | 
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4. Discussion 

Let us discuss the obtained results, giving special heed to the possible observable 
effects. 

First of all, according to (fl7|) . particle concentration near the black hole (r ~ r g ) 
is proportional to n^/voo (we should remind that and Voo are concentration 
and velocity of DMPs undisturbed by the gravitational field) . Since the annihilation 
signal intensity is proportional to the squared concentration, the signal grows by 
the factor of the order of (c/uoo) 2 with respect to the undisturbed dark matter field. 
As it has been already mentioned, the dark matter is cold; the speed of its thermal 
motion is negligible. The DMP velocity in the Galaxy is mainly determined by 
the gravitational potential of the sighting point. In the vicinity of the Sun system 
— 300 km/s. Near the Galaxy centre Woo is bigger but of the same order. 
Consequently, one can expect the signal augmentation by six orders only due to 
the dark matter concentration growth. Unfortunately, the effective region volume 
is not so large (of an order of 7n"g), but for the supermassive central black hole in 
our Galaxy r g ~ 10 7 km. 

Secondly, though the DMP speed near the black hole tends to the speed of light, 
the energy of their collisions does not excel a finite limit. Indeed, using (|6l7p one 
can make sure that the sum energy of two colliding DMPs in the center-of-mass 
system is equal: 

E = 2m ^ - J^gfg = 2m ^ =S I (21) 

where 2(3 is the angle between the velocities of the particles. The last radius, where 
the head-on collisions (with j3 = 7r/2) are still possible, is r = 2. The corresponding 
collision energy is 2V2TO. When particles get closer to the horizon, their speeds 
remain to grow, but possible (3 rapidly decreases (since the velocities of all the 
falling particles become almost parallel). When r < 2, it is determined by (|16p . 
Substituting (fT6|) to (f2"Tj) , we obtain: 

E max = 2mVr F +4 

lim E max = 2 v / 5 to (22) 

r— >1 

Thus, the maximal collision energy nowhere surpasses 2\/5 to. Practical conclusion 
from this fact is that if the threshold energy of some DMP interaction process 
excels 2-\/5 to, the process is impossible. Nevertheless, for 40 GeV particles the 
limit is - 180 GeV. 

So the energy of the DMP collisions can be very high, which can significantly 
affect the cross-sections of various annihilation channels and even change the main 
one. In particular, strong resonances at the energies around Zq and Higgs bosons 
masses can appear. Of course, it could change the annihilation products and strongly 
influence the annihilation rate. 
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Thirdly, let us consider the resulting spectrum of the annihilation signal. Un- 
der usual situations (far from compact objects) the considered way of annihilation 
(XX ~~ > 27) should produce a very narrow line: the width of photon spectral distri- 
bution is determined by the particle energy spread: 

AE ~ m DM v 2 (23) 

In the considered case there are at least three mechanisms of the line broadening and 
shifting. First, the annihilating particle collision energy significantly depends on r 
and the angle between particle velocities (see equation |2"T|) . The energy of produced 
photons varies correspondingly. Second, the centre of mass of the colliding particles 
is generally speaking not at rest. Consequently, Lorentz shifting appears. Third, 
gravitational redshift quite strongly decreases the energy of the photons generated 
near r g . We should mention that the effects of the first and third mechanisms are 
opposite in sign: so they can compensate each other. Nevertheless, one could expect 
that the annihilation line was arbitrarily wide. 

In Fig. Q] we can see that various broadening effects well compensate each other, 
and the line remains respectively narrow and unshifted. This fact is of big practi- 
cal significance. Indeed, different ways of dark matter annihilation signal detection 
from the central Galactic source (where a supermassive black hole is situated, as we 
have already mentioned) are discussed in the literature^. It is suggested to observe 
immediate annihilation products as well as the synchrotron emission radiated by 
them. However, a serious problem appears here. The concentration of not only dark, 
but of normal matter too is high in the centre of the Galaxy. The normal matter 
accretion on the black hole is accompanied by complex magnetohydrodynamic ef- 
fects which can generate high-energy particles, synchrotron emission etc. of various 
spectra and characteristics. Since the cross-section of interaction of baryonic matter 
is incommensurably greater, effectiveness of these processes is very big. Therefore, 
even if some interesting signal from the Galactic centre has been observed, one 
should prove that it is the dark matter annihilation that gives it. The case when 
the emission spectrum is wide it's a very difficult task. Whatever is the signal, to fit 
it by the ordinary matter emission is usually possible. In the case of a narrow line 
the situation is much better. It is more complicated to offer an alternative expla- 
nation to a narrow line on the energy > 10 GeV. Moreover, the energy of the line 
would immediately give us the DMP mass. Therefore, the study of the narrow anni- 
hilation line generation opportunity is a very promising problem, even if the relative 
cross-section of the channel is small. As we have seen, the DMP density near a black 
hole strongly increases. Besides, the collision energy growth can strongly elevate the 
two-photon annihilation channel efficiency. Thus, we can expect an intensive line 
formation about a black hole. A misgiving appears, however, that the line can be 
transformed into a wide band (AE/E ~ 1) by the gravitational redshift, for in- 
stance. Figure [1] shows that it does not happen. The line remains narrow enough 
(AE/ E ~ 0.05). On the other hand, the line turns out to be distinctly wider then if 
it appears far from compact objects. In the last case, the width can be estimated via 
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(|23|) which gives us AE/E ~ (woo/c) 2 ~ 10 -6 . So the line generated near a black 
hole can be easily distinguished from the line appearing under usual situations. 
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